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The fine-structure constant and the electron mass in string theory are determined by the values of 
scalar fields called moduli. If the dark matter takes on the form of such a light modulus, it oscillates 
with a frequency equal to its mass and an amplitude determined by the local dark-matter density. 

This translates into an oscillation of the size of a solid that can be observed by resonant-mass 
antennas. Existing and planned experiments, combined with a dedicated resonant-mass detector 
proposed in this Letter, can probe dark-matter moduli with frequencies between 1 kHz and 1 GHz, 
with much better sensitivity than searches for fifth forces. 


Introduction. — In string theory, the values of the fun¬ 
damental parameters, such as the fine-structure constant 
or the electron Yukawa coupling, are functions of scalar 
fields called moduli. In a typical vacuum, there are sev¬ 
eral moduli that describe geometric properties of the ex¬ 
tra dimensions of space, such as their size. The masses of 
the moduli are model dependent. Several moduli often 
remain massless as long as supersymmetry is unbroken. 

For a supersymmetry-breaking scale near a TeV, the 
moduli can acquire a mass as large as 0.1 meV or a fre¬ 
quency of 20 GHz. These rough estimates are often cor¬ 
rected by small coefficients, such as loop, logarithmic, 
and large-volume factors, that make the moduli masses 
significantly lighter mi- Moduli associated with small 
numbers, such as the electron Yukawa coupling, are also 
naturally much lighter. One may even speculate that 
the ultrasmall cosmological constant is associated with 
an ultralight dilaton whose mass is of order the Hubble 
scale [5]. 

In the absence of a general theoretical mass range for 
moduli, we will only concern ourselves with experimental 
constraints. These scalars are an excellent dark matter 
(DM) candidate when produced through the misalign¬ 
ment mechanism. In order for a scalar to be a good DM 
candidate that gravitationally clumps at galactic scales, 
it has to be heavier than 10“^^ eV [ 6 ]. For scalar DM to 
be well characterized as a scalar field, instead of individ¬ 
ual particles, its mass must be lighter than about 0.1 eV. 
Such a DM candidate—denoted by the field —can cause 
fundamental constants to oscillate in time [ 7 ]. We con¬ 
sider couplings to the electron e and the electromagnetic 
field strength 
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where Gn is Newton’s constant {h = c = ks = I 
throughout). We can identify (f with an electron Yukawa 
(electric charge) modulus if dme 7 ^ 0 (de 7 ^ 0). If 0 
constitutes the local DM energy density Pdm, it can be 


approximated by 

0(t, x) ^ cos [m 0 (t - V • X + ...)] , (2) 

rUfj) 

where |v| is the relative velocity of the DM with respect 
to Earth, roughly equal to the virial velocity in our 
Galaxy. The field oscillation occurs at an angular fre¬ 
quency equal to the DM mass, m^, and exhibits high 
fractional temporal and spatial coherence of v~-^ ^ 10 ^ 
and v~-^ ^ 10 ^, respectively, due to a low velocity disper¬ 
sion of the DM [ 8 ]. In such a background, the electron 
mass rUe and the fine-structure constant a can fluctuate 
along with (f according to 
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The size of any atom is of order l/anie, and will oscil¬ 
late if rUe or a fluctuate. For a single atom, this is a 
tiny effect. However, it is enhanced when the atoms are 
stacked, as in a solid. 

In this Letter, we show how this amplification by the 
number of atoms, in combination with resonant effects, 
can be exploited to search for scalar DM with existing 
technology already used to search for gravitational wave 
(GW) radiation. In what follows, we explain how the 
signal arises, describe the reach of existing experiments, 
and discuss future directions, including a new experimen¬ 
tal proposal. Finally, we compare with other constraints 
on scalar DM and find an improvement in sensitivity by 
several orders of magnitude. 

Experimental signature. — The effect of a DM modulus 
on a low-loss massive antenna can be captured by consid¬ 
ering the response of a harmonic oscillator. An otherwise 
free mass M on a physical, dissipative spring of equilib¬ 
rium length L, resonant frequency uj and a quality factor 
Q obeys 
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FIG. 1. Scalar field parameter space, with mass and corresponding DM oscillation frequency /(^ = m(/,/27r on the bottom 
and top horizontal axes, and couplings of both an electron mass modulus (di = dme) cind electromagnetic gauge modulus 
{di = de) on the vertical axis. Natural parameter space for a 10 TeV cutoff is depicted in green, while the other regions and 
dashed curves represent 95% CL limits from fifth-force tests (“5F”, gray), equivalence-principle tests (“EP”, orange), atomic 
spectroscopy in dysprosium (“Dy”, purple), and low-frequency terrestrial seismology (“Earth”, black). The blue curve shows 
the projected SNR = 1 reach of a proposed resonant-mass detector—a copper-silicon (Cu-Si) sphere 30 cm in radius—after 1.6 y 
of integration time, while the red curve shows the reach for the current AURIGA detector with 8 y of recasted data. Rough 
estimates of the 1-y reach of a proposed DUAL detector (pink) and several harmonics of two piezoelectric quartz resonators 
(gold points) are also shown. 


with thermal Brownian noise forces Fth in the spring and 
external noise forces Fext- In the presence of a modulus, 
the equilibrium size of the spring is oscillating in time L 
Lq cos(m 0 t). Once we define the “displacement distance” 
D = X — the influence of the modulus is revealed as a 
new force, 


M 


D+‘^D + uj^D 


—ML + Fth + Fext, (6) 


up to 0(l/(3)-suppressed force terms. The modulus- 
induced force is analogous to the tidal force caused by 
a GW [9], except that the modulus induces a monopole 
strain instead of a quadrupole strain pattern. This intu¬ 
ition can be extended to continuous acoustic systems by 
describing the modulus as a scalar GW with an effective 
isotropic Riemann curvature tensor 

( 7 ) 

where the effective strain h = —5a/a — drrielme = 
— {dfUe + ^e) inherits the coherent properties 

of the DM field oscillation as described below Eq. §. 
The response of a resonant-mass detector to modulus DM 
may thus be extracted from well-known strategies for de¬ 
tecting monochromatic gravitational-wave radiation. 

A resonant-mass detector is acoustically equivalent to 
a combination of independent harmonic oscillators, since 
the displacement from equilibrium in an elastic solid 


can be decomposed into normal modes as D (x, t) = 
y;„D„(t)u„(x). In a spherical geometry, we can take 
Un(x) = run{r), since only spherically symmetric (/ = 0 ) 
modes are excited by a scalar strain. For a sphere of 
radius R with uniform density p and longitudinal (trans¬ 
verse) sound speed q (c^), these mode functions can be 
found in Ref. m- They have resonant angular frequen¬ 
cies ujn = cikn with kn — Ti'K/R. Wc choosc a normaliza¬ 
tion such that Un{R) = 1 , and define the effective mode 
mass Mn through d^x pUn'Un' = dnn'Mn. With these 
conventions, Dn is the absolute displacement of the sur¬ 
face from the equilibrium radius Feq, and satisfies Eq. ^ 
with an effective modulus force 

Fniod,n = '^iOjO I ^ X pU^X^ — flAdQRJji (8) 
Jv 

with Mo the mass of the sphere and a coupling factor 
dr r^Un that decouples for the higher har¬ 
monics as Jn ^ n~‘^. Modulus DM can be detected if the 
force in Eq. ^ exceeds the noise forces Fth and Fext- 
Existing resonant-mass detectors .— The response of 
resonant-mass detectors to gravitational waves was first 
described by Weber mun]. Resonant-mass GW detec¬ 
tors have made great strides in sensitivity since the first 
“Weber bars” (see Ref. [13] for a historical review), so 
far culminating in a network of third-generation experi¬ 
ments mills] consisting of cryogenic, ton-scale, cylindri¬ 
cal antennas operating at around 900 Hz. Despite qual- 
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FIG. 2. Schematic of the proposed setup: a Cu-Si sphere 
whose surface displacement is monitored by a Fabry-Perot 
(“FP”) interferometer. Elements encircled by the dotted lines 
are independently suspended and isolated from vibrations. 


ity factors in excess of a million, these detectors achieve 
a sizable fractional bandwidth of 0(10%) by amplify¬ 
ing the surface displacement of the main antenna with 
a series of smaller mechanical and electrical resonators 
tuned to the frequency of the lowest longitudinal har¬ 
monic of the cylinder. The AURIGA Collaboration has 
achieved the widest bandwidth, operating at a noise level 
for 850 Hz < / < 960 Hz for an 
optimally polarized gravitational strain hij [16]. Recast¬ 
ing as a projected reach for modulus couplings at unit 
signal-to-noise ratio (SNR = 1) with 8 years of data on 
tape yields the red curve in Fig. 

Astrophysical objects provide for naturally occurring 
resonant-mass antennas. The fundamental breathing 
mode of Earth, which has a 20.46-min period and 
Q ~ 7500, was studied with a dedicated seismometer 
in Ref. [T7] over a 7-month period. The observed vertical 
acceleration noise spectrum of 7.6 x 10“^ m s“^ Hz“^/^ 
corresponds to a spherical strain sensitivity « 8.4 x 
10“^^ Hz“^/^ and a constraint \h\ < 4.5 x 10“^^. Inter¬ 
preted as a 95% CL limit on modulus DM couplings, this 
yields \de -h ^ 2.5 x 10“^ in a 1.3 x 10“^ fractional 
bandwidth around ~ 8.1 x 10“"^ Hz (black line in 
Fig. 0. Modern examinations of Earth’s higher harmon¬ 
ics [18] and crust excitations [THHH] may also be inter¬ 
esting, but likely have worse strain sensitivity. Lunar [22] 
as well as asteroseismic observations [23H25] have shown 
more promise towards detecting (quadrupole) metric 
variations; monopole strain excitations of these systems 
merit further investigation. 

Future resonant-mass deteetors. — Given the unknown 
modulus mass, there is a clear need for a wideband de¬ 
tector, or a narrow-band one with scanning ability. Pro¬ 
posals in the former category, such as “xylophone” arrays 
[26] [27] and “DUAL” detectors [28H3Q] . have been pro¬ 
posed for GW searches in the 1-10 kHz band, and would 
have excellent reach for moduli as well. The wider-band 
DUAL proposal of Ref. [30] with a capacitive read-out 
may reach sensitivity to quadrupole strains at the level 
of 10“^^ Hz“^/^. The same noise spectrum for monopole 
strains would yield a reach as represented by the pink 
curve in Fig. after 1 y of data. (This is a useful 



FIG. 3. Strain reach \h\ at SNR = 1 as a function of fre¬ 
quency / after an integration time tint = 5 x 10^ s (thick blue 
curve), consisting of a 5% fractional frequency scan by vary¬ 
ing the temperature between 4 K and 100 K in increments, 
each tshot = 10^ s long. Reach for one “shot” at 4 K is shown 
by the thin blue curve (up to the third harmonic, for clarity). 
Equivalence-principle and fifth-force exclusions are shown in 
orange and gray, while strains below the green line are natural 
for an electron Yukawa modulus with a 10 TeV cutoff. 

proxy, because although a modulus would excite higher- 
frequency modes in that setup, we expect a comparable 
sensitivity with only minor read-out changes.) 

In this Letter, we propose a scanning experiment be¬ 
cause of its simplicity and feasibility with current tech¬ 
nology. The basic experimental concept, illustrated in 
Fig-i is that of a freely suspended sphere, whose acous¬ 
tic modes can be frequency-shifted by dialing the tem¬ 
perature and detected via an interferometric read-out of 
the sphere’s surface displacement. 

By exploiting the temperature dependence of elastic 
properties, adjustment of each mode’s angular frequency 
ujn = cikn becomes possible. We propose a spheri¬ 
cal antenna of radius R = 0.3 m made out of a ma¬ 
terial such as C65500 copper-silicon alloy (high-silicon 
bronze “A”), for which q varies by about 5% below 100 
K while maintaining a high quality factor {Q ^ 10^) 
[31132]. Copper-based alloys are already used in spher¬ 
ical resonant-mass GW detectors [33l [34] for their high 
density p ^ 8 x 10^ kg/m^, sound speed q « 4 x 10^ m/s, 
and thermal conductivity. 

Brownian noise forces are broadband. By the 
fluctuation-dissipation theorem [35], their single-sided 
noise spectral density is S^pp ^ ATM^uJu/Qn around 
each mode, for a temperature T ^ Using 

Eq. (§, this translates to a near-resonance strain spec¬ 
trum 5'/^ M^Q^cf k 3 R 3 j 2 or an amplitude of 1.5 x 

10“^^ Hz“^/^(M^/Mo/i:^R^ at 4 K, with the latter 

factor scaling as 

A Fabry-Perot interferometer, schematically drawn in 
Fig. HI can measure the sphere’s total surface displace¬ 
ment X = D -h SReq through changes in the cavity length 
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^cav = 1 mm, which cause laser light to fluctuate in in¬ 
tensity on a photodiode. The shot-noise-limited displace¬ 
ment spectral density is q [l + (cj/^lcav)^], 

where the cavity bandwidth is ^^cav = 7r/2J^I/cav the 
cavity finesse is taken to be ~ 3 x 10^. We assume a 
baseline sensitivity of q « 10 “^^ m^ Hz“^ ^ X/J^‘^P, 
achievable with a laser power P ^ 1 mW and opti¬ 
cal wavelength A [36]. Converting to a strain spec¬ 
tral density yields a near-resonance strain sensitivity of 
~ 10“"® Hz-1/2(M„/MoJ„) for LOn < ficav, 
scaling as oc (above the cavity bandwidth, like n^). 
Thermal noise sources in the interferometer can be kept 
subdominant above 1 kHz with fused silica mirror sub¬ 
strates, silica/tantala coatings, and a laser beam width 
of 1 mm [371 Eg. Vibration isolation schemes with 
—200 dB attenuation exist for similar geometries [39] . 
and should be able to reduce typical seismic noise spectra 
< 10“^^ m [40] to negligible levels above 100 Hz. 

Monopole modes may be discriminated from multipole 
and other spurious modes via calibration hammer tech¬ 
niques [33] or multiple cavities [41]. 

The reach of our proposal is shown in Figs.[2&[^ The 
resonant frequencies are adjusted in fractional increments 
of 10“^ (the fractional signal bandwidth) by varying T 
in 2 mK steps. If each shot lasts a time tghot = 10^ s, a 
strain /ishot(w,Tj) ~ 

can be detected at SNR = 1 , depicted by the thin blue 
curve in Fig. [^ for a temperature = 4 K. After an 
integration time tint = 5 x 10 ^ s, a set of 5 x 10 ^ shots 

has a strain reach hint(co’) ^ ^shot(<^: less 

than 10 “^^ in a 5% band around each I = 0 harmonic 
up to n ^ 100. Better sensitivity may eventually be 
attained at ultracryogenic temperatures. Since this pre¬ 
cludes scanning of the resonant frequencies, traditional 
capacitive transducer read-outs as in Refs. Eg El or 
more advanced interferometer schemes [42] would have 
to be employed. 

Yet higher frequencies may be explored by microme¬ 
chanical resonators, for which the unfavorable scaling 
of Brownian noise with size may be mitigated by us¬ 
ing clever geometries and extremely low-loss materials. 
The proposal of Ref. [43] to detect high-frequency GWs 
in curved-plate quartz crystals is also sensitive to mod¬ 
ulus DM. An isotropic strain excites longitudinal acous¬ 
tic modes, which, due to the piezoelectric nature of the 
quartz crystal, may be picked up by an electronic circuit. 
On-resonance strain sensitivities down to 10 “^^ Hz“^/^ 
are expected for up to a hundred modes per crystal [43] . 
(Methods to increase the bandwidth are under develop¬ 
ment [44].) The modulus coupling reach of harmonics in 
two 20 mK sensors from Ref. [43] is illustrated by the 
gold points in Fig. [^ 

Non-acoustic experiments .— The phenomenology of 
light scalars with modulus couplings as in Eq. extends 
beyond the acoustic signature described above. 


Through scalar exchange, two macroscopic bodies with 
mass Ml and M 2 experience a Yukawa force with a 
range set by Its strength relative to gravity is 

<^mod - {diQ^){d 2 Q 2 ) with diQi = dm.Qme PdeQe- We 
follow the notation of Ref. [45], in which Qme {Qe) is 
the fractional amount of electron-mass (electromagnetic) 
energy relative to the rest-mass energy M of the object. 

Searches for fifth forces at a length scale L are sensi¬ 
tive to moduli with mass ^\/L. A number of exper¬ 
iments have set constraints on composition-independent 
|<^mod| < 10“^-^ for 10“^^ dSf < < 10“^-^ eV [l 6 ]. 

Rescaling by typical modulus charges Qme ^ 1/4000 and 
Qe ^ 1/500 yields approximate constraints on \dmj and 
I del shown in gray in Fig. 

The modulus force also violates the equivalence princi¬ 
ple: two test masses Mi and M 2 experience a different ac¬ 
celeration in the presence of a third one M 3 even though 
Ml = M 2 , provided that Eot-Wash 

experiment m has measured the fractional differential 
acceleration of beryllium and titanium in the Earth’s 
gravitational field to be (uBe ~ <^Ti)/<^ ~ (0.3 ± 1.8) x 
10-13. Using Qf (Qp - QP) = {-2.42 x 10“^ -3.00 x 
10-^} for i = {me,e}, one arrives at the 95% CL up¬ 
per limits \dmj ^ 1.05 x 10“^ and |de| < 2.98 x 10“^ 
for nicf) <C 1/R©, shown in orange in Eig. Bounds for 
> 1 /R© are estimated by rescaling the limit accord¬ 
ing to Ref. [48]. Lunar laser ranging sets less stringent 
constraints [45] [49]. 

Besides mediating EP-violating forces, a modulus field 
sourced by a massive body could also slightly alter fun¬ 
damental constants around it. When sourced by the 
Sun, the modulus appears as an annual modulation of 
the fine-structure constant or the electron mass with a 
known phase, and with amplitude proportional to the 
Sun’s modulus charges Qf and the ±1.65 x lO-n* an¬ 
nual variation of the gravitational potential of the Sun 
on the Earth’s orbit. The absence of such modulation 
in spectroscopy data of two different dysprosium iso¬ 
topes [50] constrains |de| to be less than 2.1 x 10“^ for 
rUff) (AU)-^. Atomic clock pairs have the capability 
to greatly improve upon this technique, and extend it to 
other couplings. 

Assuming that the field 0 comprises part or all of the 
DM density, the modulus can be probed by spectroscopic 
searches for time-varying fundamental constants [7] . Re¬ 
cently, Eourier analysis of transition energies in two dys¬ 
prosium isotopes has set the tightest constraints on |de| 
for iTiff) ^ 3 X 10-^3 eV [51], as indicated by the pur¬ 
ple curve in Eig. In the background field of Eq. Q, 
a mass M also experiences a (gradient-suppressed) force 
that may be observable in differential accelerometers such 
as free-mass GW detectors [7]. 

Black hole superradiance [52] excludes scalars with 6 x 
10-^3 eV < ^ 2 X 10“^^ eV regardless of abundance 

and the couplings in Eq. 0. unless (p has sufficiently 
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strong self-interactions [53]. Stellar cooling bounds are 
not competitive with force tests for the masses under 
consideration [54]. A detailed summary of astrophysical 
constraints can be found in Ref. [7]. 

Discussion. — The mature technology of resonant-mass 
detectors provides a new way to probe scalar DM that 
couples to the electron mass and electric charge, with a 
sensitivity beyond that of EP and fifth-force tests, over 
a wide range of frequencies. 

Resonant-mass detectors are better suited for scalar 
DM searches than free-mass interferometers such as 
LIGO, which have much reduced sensitivity to scalar 
GWs [55] because of laser phase noise. By using equal- 
length interferometer arms, this noise can be canceled 
while leaving a quadrupole GW signal unaffected. How¬ 
ever, this strategy also suppresses the scalar DM signal 
up to small gradient effects [7]. 

A scalar DM candidate frequently appears in string 
theory as a modulus, the dilaton or axions, even though 
for the latter the expected modulus couplings are far be¬ 
low our sensitivity levels [7] . Designing setups with wide¬ 
band sensitivity is crucial given the unknown masses of 
these DM candidates. Theoretically, large couplings im¬ 
ply large radiative corrections to the mass of scalar par¬ 
ticles, and thus bias towards a model-independent min- 
imum mass: [J. 

This naturalness criterion is satisfied inside the green 
bands of Fig. for a hard cutoff of A « 10 TeV. 

The technology behind resonant-mass detectors has 
steadily improved over the past five decades, and offers 
a unique opportunity to search for a well-motivated DM 
candidate over a wide range of masses and couplings in 
the immediate future. 
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